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Abstract. We overview the recently proposed mode-sum regularization prescription 
(MSRP) for the calculation of the local radiation-reaction forces, which are crucial 
for the orbital evolution of binaries. We then describe some new results which were 
obtained using MSRP, and discuss their importance for gravitational-wave astronomy. 
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The problem of including the radiation-reaction (RR) forces in the orbital evolu- 
CN ' tion of a binary is a long-standing open problem. This problem is as yet unresolved 

even in the extreme mass-ratio limit, with the particle orbiting a non-rotating black 
00 ! hole, although there has been a remarkable progress obtained from various direc- 

O ■ tions [1]. The conventional approach is to consider the fields in the far zone, and 

then use a balance argument to relate the far-zone fields to the local properties of 
G^ ! the particle. The generic failure of such approaches [2] prompted the idea to calcu- 

late the local forces acting on the particle, including the RR forces. In the following 
we discuss the RR forces acting on a scalar point-like charge, but for electric or 
^ gravitational charges the basic ideas are similar. The RR force ^^F^* which acts 

bj;)! on a point-like scalar charge q is given by [3] 

: ^^F^{t) = \\ (iz^ - u^ii^u") + + r V^Gr dr' 

r> , 16 b J-oo 

where TZ^ = R^^u" + R^aU^W^u^ — ^R^^"^^ , Rua is the Ricci tensor, is the charge's 
4-velocity, a dot denotes (covariant) derivative with respect to proper time r, 
denotes covariant differentiation, and Gr is the retarded Green's function. The first 
term is a local Abraham-Lorentz-Dirac type damping force, the second is a local 
force, which couples to Ricci-curvature and preserves conformal invariance, and 
the third is the so-called "tail" term, which arises from the failure of the Huygens 
principle in curved spacetime. The greatest problems in the calculation of the RR 
forces lurk in the tail term, because it requires the knowledge of Gr along the entire 
past world line of the charge. In addition, the self field of any particle diverges at 
the position of the particle, and the calculation of the RR forces will have to handle 
the infinities connected with the self field by providing a regularization prescription. 

Recently, Ori proposed to approach the RR problem via mode decomposition 
[4]. Ori observed, that the individual Fourier-harmonic modes of the self field are 



bounded, also for a point-like particle, although the sum over all modes diverges. 
This observation is very useful, because the calculation of the individual modes is 

relatively easy This still leaves the second, harder problem of having a rcgular- 
ization prescription to handle the mode sum. Very recently Ori suggested MSRP 
[5], which is very successful for the few simple cases to which it has already been 
applied. In what follows, we overview MSRP very briefly, and describe some of the 
recent results which were obtained using it. 

The tail part of the RR force can be decomposed into stationary Teukolsky 
modes, and then summed over the frequencies uj and the azymuthal numbers m. 
This force equals then the hmit e ^ 0~ of the sum over all I modes, of the difference 
between the force sourced by the entire world line (the bare force ^^^"^pj^ and the 
force sourced by the half-infinite world line to the future of e, where the particle has 
proper time r = 0, and e is an event along the past (r < 0) world line. Next, we 
seek a function /i^ which is independent of e, such that the series Y^i^^'^^^Fl ~ 
converges. Once such a function is found, the regularized self force is then given 
by t^'^F^ = E^C"''*'^;! - - c^M' where is a finite valued function. MSRP [5] 
then shows, from a local integration of that h^^ = a^i + + c^i~^. MSRP 
also provides an algorithm for the calculation of the functions a^jb^, and 
analytically. It has been conjectured, that for all orbits a^^ = = c^. (It was found 
to be true for all the special cases calculated so far.) If this is indeed the case, the 
tail force is given by 

t-'F^ = EC'"^.' - ^m) - d,. (2) 

Note that b^ is just the limit bare^^^oo^ ^^^^^ harepe computed using 

the Teukolsky formalism. Alternatively, b^ can also be calculated analytically us- 
ing MSRP. The only remaining problem then, is to calculate d^. Even though 
there is an algorithmic way to calculate d/^, this calculation is by no means easy. 
Although MSRP has been developed as yet only for very simplified cases, the ap- 
proach is likely to be susceptible of generalization also for more realistic cases. If 
robust, MSRP can be of the greatest importance for the calculation of templates 
for gravitational-wave detection. 

Table 1 displays the values of the MSRP parameters for the cases which have 
already been calculated. Note that the motion is not necessarily geodesic. The 
data in Table 1 may suggest the conjecture that d^ exactly equals the sum of the 
two local terms of Eq. (1) (or its analog for other charge types). If this hypothesis 
is proved to hold in general, then a truely remarkable thing happens: the full RR 
force can be calculated directly from Eq. (2) when c?^ is ignored. It should be 
emphasized that presently the support for this far-reaching hypothesis is only the 
special cases listed in Table 1. 

Next, we present some results which were obtained by application of this new 
approach. For the case of a static, minimally-coupled, massless scalar charge in 
Schwarzschild, the self force is known to equal zero [11]. For a static electric charge 
q in Schwarzschild the self force is known to be purely radial and to be given by 



TABLE 1. Values of the MSRP parameters. All spacetimes are spherically symmetric. For all 
the cases = = c^. The charge's spin is s. A * denotes the cases for which and were 

inferred indirectly, and a f denotes cases where the values of &^ were corroborated numerically. 
Here, = 2-Fi(l/2, 1/2; 1; w^) and h = 2^1(1/2, 3/2; 1; i;^), = -(d^/ dtfr^ / gtt, is the 
4-acceleration, o? = iiau", and fc^ = ^(it^i — u^a?). 
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fr = g2Mr-3(l - 2M/r)-^/2 [12]. These results were recovered using MSRP in 
[7]. Note that for these two simple static cases the solution for the modes can be 
obtained analytically. In general, however, this is not expected to be possible, and 
the solution can be obtained only numerically. 

The case of a scalar charge q in uniform circular orbit around a Schwarzschild 
black hole was recently considered in [9] . The RR force was calculated numerically 
without any simplifying assumptions, such as far field or slow motion, and the 
solution is fully relativistic. Both the temporal and the azimuthal, dissipative 
components and the radial, conservative component of the RR force were computed. 
Figure 1 displays the behavior of the radial component of the RR force for both 
geodesic and non-geodesic orbits. In the slow motion and far field limits the force 
is repulsive, and behaves like ^^Fr ~ g^M^fi^/r^. However, in strong fields the 
force grows faster, and for fast motion is changed from repulsive to attractive. 
This expression for the radial, conservative RR force may be very important for 
the detection of gravitational waves, and also for gravitational-waves astronomy. 
The conservative radial force causes an additional precession of the pcriastron of 
the particle's orbit, and thus induces a change in the frequency and phase of the 
emitted radiation [13]. Although the radial self force has been obtained only for the 
simple case of a point-like scalar charge, the result indicates that one can expect a 
non-zero periastron precession also for a small mass. However, the magnitude of the 
effect will very reasonably depend on the type of the charge. A large-magnitude 
effect can cause the entire search algorithm to fail in the very detection of the 
signal (depending also on the size of the template library), and a small effect will 
introduce errors in the parameters of the observed binary, namely, the wave form 
would fit the template of a system with parameters different from the parameter 
of the actual binary. Note that the conservative force depends not only on the 
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FIGURE 1. The radial RR force acting on a scalar charge in uniform circular orbit around a 
Schwarzschild black hole. Top panel (A): Fr as a function of r/M for geodesic orbits. Bottom 
panel (B): Fr as a function of the angular velocity f2 in units of the Keplerian angular velocity 
CIk, when the orbit is at r = 4M. 



radiative modes of the field, but also on the non-radiative modes [13]. 
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